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Abstract 
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I Information metrics give lower bounds for the estimation of pa- 

■ rameters. The Cencov-Morozova-Petz Theorem classifies the mono- 
. tone quantum Fisher metrics. The optimum bound for the quantum 
I estimation problem is offered by the metric which is obtained from the 

■ symmetric logarithmic derivative. To get a better bound, it means to 
CN ■ go outside this family of metrics, and thus inevitably, to relax some 

I general conditions. In the paper we defined logarithmic derivatives 

^ ■ through a phase-space correspondence. This introduces a function 

^ . which quantifies the deviation from the symmetric derivative. Using 

^ I this function we have proved that there exist POVMs for which the 

■ - - ■ new metric gives a higher bound from that of the symmetric deriva- 
tive. The analysis was performed for the one qubit case. 
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1 Introduction 



In the classical parameter estimation theory a fundamental role is played by 
the Cramer-Rao inequality. According to the latter the variance-covariance 
of an unbiased estimator is bounded from below by the inverse of the Clas- 
sical Fisher Information Matrix, which is constructed from the logarithmic 
derivatives of the parametric statistical model. The Fisher Information Ma- 
trix provides a metric on the information manifold of the probability distribu- 
tions of the model. A basic property of such a metric on a statistical manifold 
is its monotonicity under stochastic (information loosing) transformations. 
This means that the distributions are closer w.r.t. such metrics or become 
less distinguishable under such transformations. The well known theorem 
of Cencov [5] states that the Classical Fisher Metric is the unique metric 
with the above monotonicity property. In the Quantum Parameter Estima- 
tion Problem there are new aspects in the parameter estimation inequalities. 
First a parameter of a quantum state is estimated using appropriate quantum 
measurements which provide classical probability distributions. These dis- 
tributions give the classical Fisher metrics with corresponding Cramer-Rao 
inequalities. Second the question of the existence of metrics which are inde- 
pendent of the measurements and characteristic only of the quantum states 
is related to the existence of Quantum Fisher Metrics with the monotonic- 
ity property under completely positive maps of the states. The Theorem of 
Cencov- Morozova-Petz [T^l [13] gives the general form of all monotone met- 
rics. With respect to these metrics there is a hierarchy of Cramer-Rao type 
inequalities which provide bounds for the variance-covariance of unbiased es- 
timators. While the highest bound is given by the Classical Fisher metric 
of the measurement probability, the quantum monotone metrics give lower 
bounds [9J. The highest among them is given by the symmetric logarithmic 
derivative (SLD) and the lowest by the right logarithmic derivative (RLD). 
There is an extensive literature on the saturation of these inequalities for fi- 
nite ensembles and in the asymptotic regime [21 [H [15] . We posed a question 
in the opposite direction. While the saturation of the inequalities cannot 
be achieved for all measurements and all types of quantum states, we asked 
whether there exist higher bounds for a class of measurements that give a less 
informative estimation. This would be useful in classifying the effectiveness 
of certain measurements which may actually be realized. Existence of higher 
bounds means less accuracy in parameter estimation for certain measure- 
ments. Given the fact that the SLD gives the highest bound which, anyway, is 
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smaller than the classical Fisher bound, a higher bound, in the spirit of quan- 
tum metrics, can be derived if one goes outside the family of metrics resulting 
from the Cencov-Morozova-Petz Theorem. Since the non-uniqueness of the 
quantum information metrics is related to the non-commutativity of the log- 
arithmic derivatives, it seemed natural to construct logarithmic derivatives 
indirectly through a phase-space correspondence. That is, for the given quan- 
tum state, we took the logarithmic derivatives of the corresponding Husimi 
function and mapped them back as operators in the original Hilbert space. 
Using these operator valued objects we derived a metric. But, as expected, 
at least for the case of one qubit, the form is not that of the Petz family of 
metrics. This means that our metric is not manifestly monotone. Neverthe- 
less we investigated what could be its role in the hierarchy of Cramer-Rao 
type inequalities. We introduced a function which in some sense quantifies 
the deviation of the new metric from the SLD metric. Using this function 
we proved that there exist measurements for which the mean square error 
is bounded from below, with a bound even higher than the classical Fisher 
bound. This result could be useful in the investigation of the effectiveness of 
certain estimation procedures. 

In the paper we first give a short introduction on the classical and quantum 
estimation problems with the corresponding inequalities. Then we present 
the necessary tools for the phase-space correspondence. Next we construct 
the new logarithmic derivatives and the new quantum metric. Using this we 
state a theorem concerning the new bound and apply it to an example of one 
qubit. The proof is given in the appendix. Finally we discuss open problems 
and further possibilities of this approach. 

2 The classical and quantum parameter esti- 
mation problem 

Here we present only the necessary concepts in order to establish the notation. 
We refer to the bibhography for the details [H [9] . Let 

S = {p^ = p{x;0\^=[e,-,C]^^} (1) 

be a parametric family of probability distributions on X. This is an n- 
dimensional parametric statistical model [1]. Given the N observations Xi, xjy, 
the Classical Estimation Problem concerns the statistical methods that may 
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be used to detect the true distribution, that is to estimate the parameters ^. 
To this purpose, an appropriate estimator is used for each parameter. These 
estimators are maps from the parameter space to the space of the random 
variables of the modeL The quahty of the estimation is measured by the 
variance -covariance matrix = [v^^] where 

vi = E^mx)-e){i^{x)-e)] (2) 

Suppose that the estimators are unbiased, namely 

e£{x)] = e G S (3) 

Then a lower bound for the estimation error is given by the Cramer-Rao 
inequality 

VS > G{0-' (4) 

where G(0 = [g^,m 

g,jiO = Edm^;Wix;0] (5) 
the Classical Fisher Matrix with 

k = K^'^ = logp{x; ^) (6) 

the score function. As it has been shown the Fisher Matrix provides a metric 
on the manifold of classical probability distributions. This metric, according 
to the theorem of Chencov is the unique metric which is monotone under 
the transformations of the statistical model. This means that if the map 
F : X y induces a model Sp = {q{y] ^)} on 3^ then 

Gf{0 < GiO (7) 

That is, the distance of the transformed distributions is smaller than that 
of the original distributions. Thus monotonicity of the metric is intuitively 
related to the fact that in general we loose distinguishability of the distribu- 
tions after applying any transformation to the information. 
The case of Quantum Estimation Problem differs in two fundamental aspects 
from the classical one. First the parameters to be estimated are classical pa- 
rameters of the apparatuses that prepare the quantum system to a given 
state. These parameters are measured through an interaction of the system 



4 



with appropriate measuring devices. The measurements give classical dis- 
tributions which must be estimated as in the classical case, and thus the 
estimation errors obey the corresponding Cramer-Rao inequalities. The ef- 
ficiency of this procedure depends on the effectiveness of the measurements 
used. The immediate question is what is the ultimate bound for the esti- 
mation error. This bound must depend only on the state of the system and 
as it turned out it is given by the inverse of a quantity termed Quantum 
Fisher Information Metric [H [9]. Let 11 = {7r{x)} be the measurements used 
to estimate the state of the system. Then the probability distribution of 
the measurements is p{x\^) = tr[p^7T{x)]. The corresponding classical Fisher 
information metric Gn(0 satisfies the inequality 

GuiO < GUO (8) 

where Gg^siO is the Quantum Fisher Information Metric of [ISIIH], while the 
index s means symmetric and indicates the way this metric is constructed, 
that is the way the logarithmic derivative is generalized to an appropriate 
operator. More specifically this operator valued quantity is defined implicitly 
by the equation 

^ = liP^L, + L^P,) (9) 
and then Gg^siO is defined by the equation 

9^,J,,A0 = Hpi,\{UL, + L,U)) (10) 

From now on we write gq^Cj as gsLoiO denoting that it is constructed 
through a symmetric logarithmic derivative. The SLD metric is monotone 
under completely positive maps, but it is not unique. This fact is the second 
fundamental difference from the classical case. There are various approaches 
for the definition and construction of Quantum Fisher Information Metrics 
[H El [121 [13] ■ The Cencov-Morozova-Petz Theorem [121 US] gives a general 
characterization of monotone metrics using operator monotone functions [121 
[T51 [TT] . For these metrics we have the inequalities 

V^e(0 > GuiO'' > GsLDiO'' > GMiO'' > GRLDiO'' (11) 
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where Gm(0 is any monotone metric and GRLoiO is the metric correspond- 
ing to the right logarithmic derivative defined as 



(12) 



For the case of one qubit that we analyze here, the monotone metrics have 
the general form [121 113] 



where g{t) = and f{t) satisfies f(t) = tf(t^^). The above inequalities 
satisfied by the quantum information metrics lead to the conclusion that, 
as far as the estimation problem is concerned, only the SLD metric has any 
practical use, though there are some cases where other metrics play a role [9]. 
A great part of the literature on the estimation problem deals with the ques- 
tion of saturation of these inequalities mainly in the asymptotic case [HISIIB]- 
In the present work we pose a question in the opposite direction. We ask 
whether there exit measurements for which there is a higher bound in the es- 
timation error. In the sequel, after introducing our phase-space construction 
of the logarithmic derivatives, we prove the existence of such measurements 
with a higher lower bound even higher than that of the classical Fisher metric 
of the corresponding measurement. 

3 Phase-space correspondence 

Since our interest is on the question of the existence of measurements with 
higher lower bounds of estimation errors, we restrict our analysis to the one 
qubit case. Then the Hilbert space is 2-dimensional and the information 
manifold 3d |3]. The state of one qubit is given by the density operator 



where a = (cti, (T2, as) and h = {sinOcoscp, sinOsincj), cos6) with < r < 1 , 
O<0<7r,O<0< 27r. Explicitly we have 




(13) 




(14) 



p(r, e, 0) 
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1 



1 + rcosO re '^'f'sinO 
re^'^sinO 1 — rcosO 



(15) 
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r,6,(f) are the coordinates of the information manifold for the qubit. The 
dynamical symmetry group for the system is SU{2) and the phase space is 
the unit sphere = SU{2)/U{1). Let n{e,(f)) be a point on S\ The Weyl 
Rule 

FA{n; s) = tr[AA{n; s)] (16) 

is a map of an operator A acting on the Hilbert space to a function on the 
phase space. {A(i7;s)} is a parametric family of operators defined on the 
phase space, and called Stratonovich-Weyl Kernels [1]. The inverse map is 

A = J dfx{n)FA{n;s)A{n;-s) (17) 

where here dn{Q) = sin6d6d(j)/2TT. We have 

1 I 
A(^,0;s) = y2^J]<l/2,l/2;/,O|l/2,l/2>-^ J2 DirnYi*^{e,(j)) (18) 

1=0 m=-l 

where 



1/2 



A™=J^^ < l/2,A;;/,m|l/2,g> |g>< A;| (19) 

k,q=-l/2 

or explicitly 

, . If l + 'i^cosO 3^e-''''sin9 , . . 

'^(^''^'^) = 2|^ 3^e^W l-3^co.^ I 

For A = p we get the quasiprobability distributions (QPD) P{fl; s) = Fp{Q; s). 
For s = — 1 we get the Husimi Q-function, which is non-negative and thus a 
probability distribution. We have 

Q{r,9,(f); 61, (pi) = tr[p(r, 6^, 0) A(6'i, 01, — 1)] = -{l+rcosdiCos6+rcos{(j)i—(j))sin9isini 

(21) 

Now having a genuine probability distribution associated directly to the state 
of the system, and not specific to a measurement, we may construct a Fisher 
metric. 

9ij,Husimi = EQ[dilogQ{^; n)djlogQ{^; n)] (22) 
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This was proposed as a way to avoid the non-commutativity difficulties of 
the direct Hilbert space constructions [SI [7] . In [7] a classical Fisher metric 
was constructed for a squeezed state. As it turned out, this metric was 
proportional to the SLD one. In [H] using this phase space approach a 
metric was given in the form 

-2r-log(^) , 2r + (1 - r^og (^) , 

This metric is manifestly non-monotone, as it differs from the standard form 
of the monotone metrics fl2[ [13] . The conclusion is that, though the con- 
struction of the classical Fisher metric through the phase-space distribution 
is very natural, it turned out that either it gives no new metric, or it gives 
one manifestly outside the monotone family. Motivated by these results we 
looked at an alternative route, that of a construction of quantum logarithmic 
derivatives through the phace-space correspondence. 



4 A new quantum information metric 



We propose the following equally natural construction: Starting from the 
quantum state, construct the phase-space distribution (Husimi function) and 
take partial derivatives with respect to the relevant parameters of the model. 
These are functions on the phase-space. Then using the Stratonovich-Weyl 
Kernels map these functions back to Hilbert space operators. Finally use 
these operators to construct, in the standard way, the Quantum Fisher Met- 
ric. We present our results for the one qubit case. 
We have for the new operator logarithmic derivative 

'"2" sin 9id(t)id9i 

(24) 



U = L,{r, e,(j))= r r dilogQir, 6, 0, 9^, 0i)A(0i, 0i, 1)^ 

J 01=0 J (1)1=0 

We get 



1 / f '^-^' W ( r-3cose -3e-''^sine > , . 

2^ V V T+^J J V -Se^'^stne r + 3cose ' ^^^> 



3 f ^ — f —sinO e ^'^cosO 



_ 2r + (1 - r'Vos — — ) (26) 
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Prom these we have 

= ^ (-^^ - "^'^i!? ^ - ^'vo, 

(28) 

This is our new quantum Fisher metric, which is manifestly non-monotone. 
Using this metric we construct a new bound for the estimation error which 
for some measurements is even higher than that of the corresponding classical 
Fisher bound. 



5 Existence of a new bound 



We introduce three operators which quantify the deviation of our new metric 
from the SLD metric 

^ ^'^^ ^^^^ 



^ ^ 4r^ + 2rlog (^) / i o \ 6r + (3 - r^)log (i^) / cos9 smOer^'f 



4r3 V 1 y 4r3 V sm^e'* -cos9 

(30) 

6r - Ar^ + 3(1 - r^)log (i^) f sind -cosde 



~ 8r2 ' (, -cosOe''^ -sine ' ^^^^ 

_ 6r - + 3(1 - T^)log f o mn^e^'^ n 



8r2 V -isinOe^'t' 

These operators are involved in the construction of the bound. Here we define 
the main quantities and present the main results. The detailed proof of the 
theorem is presented in the appendix. 
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5.1 The General Bound 



We suppose that we approximately know the value of the parameter ^, with 
an error e and that, to be general, the POVM measurement depends on 
i.e. n = n(^, C,'), where ^' = .^ + e. Then the classical probability distribution 
to be estimated is qe{C,\0 = t^lpiO^i^^O] We assume that the estimation 
is unbiased, namely 



Let = (2/1,1/252/3) and Z* = (2:1, 2:2, 2:3). We define the quantities 



(33) 



T^ = iJ2 vM' - e) n = m) = E 'M^^ 

i=i j=i 



(34) 



^i(e) = e /tr(p(0^^^^|^Me 



We have 

Theorem 1. Let Vg^ = {v'l} 



the covariance matrix w.r.t. qe{^\C)- 
Then we have 



(35) 
(36) 



(37) 
(38) 



[Y%J) (Z'GZ) > 



Y'{A + B)Z~tr T^nm. i + ^)di 



(39) 



The proof can be found in the Appendix. It is based on a similar Theorem 
of [TO]. 

The above inequality is general, and it holds for an arbitrary measurement. 
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The crucial point is that it depends on the operators h which depend only 
on the state, and not on the measurement. This means that for a given state 
we may extremize w.r.t. the measurements in order to uncover a possible 
new lower bound for a class of measurements. Since the general case is quite 
involved, and since our purpose is to show that there exist such bounds, we 
analyze in the next section the case of one unknown parameter. 



5.2 One Dimensional case 

For the one qubit case the parameters are ^ = (r, 9, 0). Considering to be 
the single unknown parameter the theorem states that 



where 



'(/>+€— TV 

2 



r-<p+e+-K 

a{(l)) + h{(j)) -tr ( - ) /i3n(0, (j) + e)d(i) 



(40) 



» = lis (2'- + (1 - '■>9T:p7j («) 

K,) ^ ,-j22iir[m?^^^)4 (43) 
A general POVM in the present case has the form 

''^'^''^^ ^ x'\^A)-iy'\l<t>) x^\lct>) ^ ^^^^ 

Then we get 

a{(f)) + b{(f)) — 2h^r{r)sin9 I icj) — {sin(j)x^'^ + cos(j)y^^)d4> 

(45) 

where 

6r - 4r3 + 3(1 - r2)% (i-^) 
hsr = g^. ' '^'^'^ (46) 

We note that the quantity inside the absolute value, according to the proof of 
the Theorem 1, is real. Now for measurements independent of the estimated 
parameter we can state the result as a Corollary 
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Corollary 2. Let the measurement be independent of the estimated param- 
eter (f), then since a(0) =0 and 6(0) = 1 the inequality simplifies to 

Vq,g - (^^~ 2h3r{r)sine J (j> - (p^ {sincpx^'^ + cos4>y^'^)d^^ (47) 

We observe that the bound concerns a family of estimators with the same 
nondiagonal elements. 



5.3 Maximizing the bound 

We want to maximize the right hand side of the inequality. This means to 
find conditions on a class of measurements that maximize the quantity 



C = ^a(0) + b{(f)) - 2h3r{r)sin9 J (j) - <f>^ {sin(f)x^'^ + cos(f)y^'^)d4> 

(48) 

Our purpose is to prove the existence of measurements for which the bound 
is higher than that of the classical Fisher metric. We impose conditions 
which the measurements must satisfy and which make them to have a very 
restricted form. This is only to prove our point, that is to prove that the 
class of measurements for which there exists a higher bound is not empty. 
A better approach would be to determine the most general form which has 
the property of being bounded with this higher bound. The non-negativity 
of the POVM gives the non-negativity of its eigenvalues, namely 

Ai = i (^X'' + X^^ - ^(a;ll _ a;22)2 + 4(^12)2 + 4(^12)2^ > Q (49) 
A2 = ^ (x^^ + X^^ + V(xll - ^22)2 + 4(^12)2 + 4(yl2)2^ > Q (50) 

The unbiasedness condition of Eq 33 with its LHS changed to -|- is 

[^^'^'' fl .\fl+rcose 11 l-rcosO 22 ■ n t ^ 12 ■ ^ i2\\ 
I 10— (0 + ejl I X H X + rstnO [coscpx — stncpy ) I 

Jip+e-n ^ / V 2 2 J 

(51) 

where we have made the following general assumptions and specific choices 
for the matrix elements of the POVM 

(0, + e) = x^^ (0, + e) > (52) 
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/•4>+e+iT 

J (p+e—TT J (p+e—TT 



/ x^\(i),(j) + e)d(f) = l 

J (j>+e—TT 

cl)+e+iT 



(53) 
(54) 



We assume further that a;^^(0, + e), as a function of is symmetric w.r.t. 
+ e. Then if we choose 

= {tan(t))y^'^ (55) 
the condition given by Eq 51 is satisfied and we also have 

sin(f)x^'^{$, + e) + cos(l)y^'^{$, + e) = -?/^^(0, + e) (56) 



cose 



Then the quantity to be bounded becomes 



+ e) rf0 



C = (^a(0) + 6(0) - 2hr{r)sine ^'^^ (0 - 0) (^l/''(0, 

(57) 

With the given assumptions and special conditions a bound for C comes from 
the nonnegativity of the eigenvalue Ai. We get 



a;^^(0,0 + e) > 



7/12(0, + e) 



coscj) 



(58) 



Assuming further that 



yl2 


>, + e) 




COS0 




;0,0 + e) 




coscj) 




>,0 + e) 




coscp 




c{<t)) 



> for < + e 
= for = + e 
< far > + e 



we have 



C < [c(0) - 2/i3,(r)sm^ {h - h)? 
13 



(59) 
(60) 
(61) 
(62) 

(63) 



where 



^1=/ ' (U-{cl> + e))x''{^,cl> + e) + e-^y''{lcj> + e)]d^ (64) 

72= / U-0)x^^(0,(/) + e)# (65) 

Thus for the selected class of measurements we have a lower bound for the 
variance of the estimator which is different from that given by the SLD. The 
usefulness of this result is that this bound, for this class of POVMs is higher 
than not only the bound provided by SLD, which is higher than all those 
provided by the monotone metrics, but it is higher even from the bound of 
the corresponding Classical Fisher metric, as it is shown in the next section. 

5.4 An application 

To make the bound explicit we choose the form of so that an analytic 
integration is possible. We take 



and from Eq 58 considered as equality, and with the appropriate signs. 
With this choice we get a(0) =0 and 



(67) 



(^) 



For the classical Fisher metric for the measurement we have 

9o,Fisher{r, 9,(j)) = J ^ J^^^^ (po(0|0; r, ^))^ Po(0|0; r, 9)d(p (68) 
where 

poim r, 0) = tr [p{r, 0, 0)n(0, 0o)] (69) 

Now in order to plot the different bounds we have chosen a value for 
so that cos4> < and 9 — '!r/2. In Figure 1 we plot B^ax — Cmax/g and the 
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Figure 1: Comparison of the different bounds for the variance of the estimator 
for (p. For the plot, was chosen so that cos4> < and 9 = tt/2, Bmax — 
Cmax/g- BsLD, Brld, Bpisher are the corresponding quantities for SLD, RLD 
and classical Fisher metric. The x-axis is the parameter r of the state which 
is related to the puriness of the qubit. 
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corresponding quantities Bsld, Brld, Bpisher for SLD, RLD and the classical 
Fisher metric versus the parameter r, which is related to the puriness of the 
state. We see that the new bound is higher than all the others. The deviation 
becomes bigger for higher values of r. 



6 Discussion 

We have posed the question whether, for certain classes of quantum estima- 
tors the lower bound for the variance of the estimated parameter is higher 
than those given by the quantum Cramer-Rao inequalities. The current 
literature deals with questions of saturation of these inequalities, that is 
which measurements are optimal with respect to their accuracy of estima- 
tion. Here we asked whether there exist "bad" measurements for which the 
error is bounded from bellow with higher bounds. This would be useful in the 
classification of measurements with respect to their effectiveness. The high- 
est possible bound given by the monotone quantum metrics is that of the 
Symmetric Logarithmic Derivative. This is majorized by the classical Fisher 
bound associated with the probability distribution of the measurement. To 
go beyond these bounds means to construct a metric outside the family 
of manifestly monotone metrics. To do that we constructed the operator 
logarithmic derivatives indirectly through the phase space correspondence. 
That is we mapped the state to the Husimi function, we took the logarith- 
mic derivative and we mapped them back to Hilbert space operators. The 
related metric is manifestly not monotone. Then using an operator which 
quantifies the difference of this derivative from the symmetric one we were 
able to derive a bound for the variance. This depends both on the state 
and on the measurement. We were not able to determine the general class 
of the measurements for which the new bound is higher than the previous 
ones, but using an explicit example we showed that this class is not empty. 
Now observing that, under the conditions we have imposed on the measure- 
ment, we have that one of its eigenvalues is zero, namely Ai = 0, we see that 
the measurement belongs to the class of sharp measurements [15]. Further 
work is needed to determine the class of all measurements which have this 
property of bounded accuracy and, most importantly, to analyze the physical 
realizability and nature of such measurements. 
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A Proof of Theorem 1 



The proof is based on similar techniques used in [lO]. First the unbiasedness 
condition is 

(70) 



Taking the derivative w.r.t. ^'^ we have 

( / H- tr (/ «.^n(« H- ) - SI (71) 

On the other hand we have 

etr I / + etr ( f ^n(« + 1 ^ (72) 



Subtracting these equations we get 

tr (I {i^ - e) ^n(^; ^ + e)di^ =Bi + Ai (73) 
or using Equ. 29 

tr (1 (e^- - e) (/ifc + ^(pLfc + L.pij n(e, e + e)t/e) =Bi + Ai (74) 
Now using a Lemma in jlOj we can write this equation as 



trl / Uj-^AhkU{^,^ + e)dn + 



Re 



Bi + Ai 



(75) 
(76) 



Multiplying the last equation with the real numbers yjZk and summing over 
j,k we get 



Re 
where 



tr / T^pTLU{^,^ + e)d^ 



Y\B + A)Z-tr / T^ThU{^,^ + e)d^ (77) 



(78) 



k=l 
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Thus 



tr / pTiJi{l^ + e)T^di 



> 



Y\B + A)Z - tr / T^nm, ^ + e)di 



and from Schwartz inequahty we have 



tr J pT^Uii^ + e)T^di^ (^tr J pnud ^ + e)ndi^ > 



Y\B + A)Z - tr / T^T^m. i + ^)di 



(79) 



(80) 



(81) 



This gives the inequality that was to be proved. 
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